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Abstract 

We illustrate the mass and charge renormalization procedures in quantum 
field theory using, as an example, a simple model of interacting electrons and 
photons. It is shown how addition of infinite renormalization counterterms to 
the Hamiltonian helps to obtain finite and accurate results for the S'-matrix. In 
order to remove the ultraviolet divergences from the Hamiltonian, we apply the 
Greenberg-Schweber “dressing transformation” and the Glazek-Wilson “similar¬ 
ity renormalization”. The resulting “dressed particle” Hamiltonian is finite in 
all orders of the perturbation theory and yields accurate S'-matrix and bound 
state energies. The bare and virtual particles are removed from the theory, and 
physical dressed particles interact via direct action-at-a-distance. 


1 Introduction 

Consistent nnification of relativity and qnantnm mechanics remains an unsolved theo¬ 
retical problem in spite of many efforts applied to its solution in the 20th century. The 
fundamental difference between relativistic and non-relativistic physics follows from 
the famous Einstein’s formula E = me?. This formula, in particular, implies that if a 
system of particles has sufficient energy E of their relative motion, then this energy 
may be converted to the mass m of newly created particles. Generally, there is no limit 
on how many particles can be created in collisions, so any realistic quantum mechanical 
description of high-energy systems should involve states with any number of particles 
from zero to inhnity. The number of particles is not conserved during time evolution. 
The most familiar example of such a behavior is the emission and absorption of light 
(photons) in electrodynamics. 
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First attempts to describe relativistic quantum systems were undertaken immedi¬ 
ately after creation of the formalism of quantum mechanics in 1920’s. A quantum 
theory of the electromagnetic field was constructed by quantization of the classical 
Maxwell electrodynamics. In lowest perturbation orders, this theory agreed well with 
experimental results. However, perturbative calculations for the F-matrix did not work 
in higher orders, in particular, due to ultraviolet divergences. 

The way to calculate the S'-matrix in QFT accurately in all orders was provided 
by the renormalized QED formulated by Feynman, Schwinger, and Tomonaga in the 
late 1940’s. However this approach created a host of other problems. According 
to the prevailing interpretation, the creation and annihilation operators present in 
the Hamiltonian and Lagrangian of QED correspond to bare particles having inhnite 
masses and charges.^ However, the bare particles have never been directly observed 
in experiments. They are believed to be surrounded by clouds of virtual photons and 
electron-positron pairs, thus forming complex objects called dressed particles. The 
dressed particles are supposed to be the eigenstates of the full Hamiltonian. They 
have finite experimentally observable masses and charges. The problem is that the bare 
particle Hamiltonian of QFT is formally infinite. Although, these infinities cancel out 
when the S'-matrix is calculated, the Hamiltonian is useless if one wants to calculate the 
time evolution or to hnd wavefunctions of bound states via diagonalization procedure. 

Two lines of research were initiated to cope with these problems. The dressed 
particle approach was suggested by Greenberg and Schweber [T]. Their goal was to get 
rid of the bare particles and to express the entire formalism of QFT through observable 
dressed particles only. The operators for dressed particles were sought as unitary 
transforms of bare particle operators. From another direction, Glazek and Wilson [2j 
introduced the similarity renormalization formalism. The idea was to apply a unitary 
transformation to the Hamiltonian in order to ensure that interaction potentials rapidly 
fall off as functions of energy differences, and to guarantee that all loop integrals are 
convergent. The similarity and dressing transformations are strikingly similar. They 
can be combined into one unitary similarity-dressing transformation which achieves 
two goals at once: the theory can be expressed in terms of real dressed particles and 
the Hamiltonian can be made hnite in all perturbation orders. At the same time, 
the accurate and well-tested S'-matrix of the renormalized theory remains intact. This 

^Infinite quantities appear frequently in the renormalization theory. In order to avoid them in 
practical calculations, one often uses regularization, e.g., cutting off all integrals at large integration 
momenta. In a regularized theory all quantities become finite. However, the regularization is an 
artificial trick, and in order to get rigorous and accurate results one should take the infinite limit 
of the cutoff momentum. As will be discussed later in the paper, the goal of renormalization is to 
introduce certain cancellations between (finite) regularized quantities, so that when the integration 
cutoff is lifted the physically relevant parameters have well defined finite limits. Without indicating 
this explicitly in calculations, we will always assume that the regularization and the taking the cutoff 
to infinity steps were properly executed. A quantity will be called finite (infinite) if it has (does not 
have) a finite value in the limit of infinite cutoff. 
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combined approach was developed in refs. wu and dubbed the relativistic quantum 
dynamics, or RQD. 

In this paper we will discuss the reasons why the renormalization difficulties (e.g., 
the absence of a well-dehned Hamiltonian and unsatisfactory treatment of the time 
evolution) persist in current relativistic quantum held theories. To avoid unnecessary 
mathematical complications, we will be working with a simple model theory which, 
nevertheless, shares some important features with QED. In section 9 we will explain 
how the similarity-dressing transformation of this model theory leads to a well dehned 
(perturbatively) hnite Hamiltonian that can be used for all kinds of quantum mechan¬ 
ical calculations (5'-matrix, bound states, time evolution, etc.) without the need for 
renormalization. In addition, the RQD formalism does not use the dubious notions of 
bare and virtual particles. 

2 Bound states, time evolution, and scattering 

In this paper we will be concerned with three types of phenomena that are normally 
studied in physical experiments: bound states, time evolution, and scattering. These 
three areas account for the most of experimental information available about fun¬ 
damental particles and their interactions. The key theoretical quantity involved in 
quantum mechanical description of these phenomena is the Hamilton operator H. The 
energies En and state vectors |T)„ of bound states can be found as eigenvalues and 
eigenvectors of the Hamiltonian 


= ( 1 ) 

The development of the state vector |\k) from time t' to time t is described by the 
time evolution operator exp{iHt) 

|^(t)) = (2) 

If the eigenvalues En and eigenvectors |\k)„ of the Hamiltonian are known and the 
initial state is represented as a sum (and/or integral) over the basis states 


n 

then the time evolution can be calculated as 


n 


(3) 

(4) 
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Unfortunately, in most cases, the full spectrum of the Hamiltonian is not known, and 
the time evolution is difficult to predict. 

This difficulty, however, is not that disappointing, because experimental observa¬ 
tions of the time evolution in subatomic world are even more difficult than calculations. 
Most experiments in high energy physics are performed by preparing free particles or 
their bound states (like hydrogen atoms or deuterons), bringing them into collision, 
and studying the properties of free particles or bound states leaving the region of col¬ 
lision. In these experiments, it is not possible to observe the time evolution during 
interaction: particle reactions occur almost instantaneously and one can only register 
the reactants and products which move freely before and after the collision. This gives 
a lucky break for theoreticians: In such situations the theory is not required to describe 
the detailed dynamics of particles during the short interval of collision. It is sufficient 
to provide a mapping of free states before interaction onto the free states after the 
interaction. To describe scattering experiments, one needs only the formula for the 
time evolution from the remote past T -C 0 to the distant future t ^ 0 


where the S-operator is defined by 

S = lim lim (6) 

t^OO t'^ — OO 

One can read the right hand side of eq. from right to left as a sequence of three 
steps: (i) the non-interaction evolution of the system from time t' in the past to 0; (ii) 
the sudden jump at t = 0 described by the 5'-operator; (iii) the free evolution from 
f = 0 to the future time t. 

The most effective technique available for calculations of the S'-operator is the 
perturbation theory which can be written in many equivalent forms.^ The Dyson 
time-ordered expansion provides the most economical expressions that can be encoded 
in familiar Feynman diagrams. However, for the discussion in this paper, we found 
more useful two other perturbative expressions which differ from the Dyson’s formula 
only by re-shuffling the terms. The “old-fashioned” formula for the S'-operator is 

/ + 00 r + OQ pt 

V{t) dt - / V{t) dt / U(F) dt' + ... (7) 

-oo J—oc J—oo 

where V is the interaction part of the total Hamiltonian H = Hq + V, 

^In this paper we are not discussing the complicated issue of the convergence of perturbative 
expansions. We will assume that all perturbative series do converge. 
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( 8 ) 


V{t) = 

and the factor in the limit e —> 0 serves for adiabatic switching the interaction 
on and off. Operators with t-dependence determined by the free Hamiltonian Hq as in 
eq. (jHI) will be called regular. Using convenient symbols for f-integrals 

= t Y{t')dt' 

J — OO 

/ + 00 

Y{t')dt' 

-OO 

formula dZD can be written compactly as 



where 


5 = 1 + S(f) 


(9) 


E(«) = iV(t) - V(t)V(t) - iV(t)V(t)V(t) + V(t)V(t)V(t)V(t) + ... (10) 

Another equivalent expression for S was suggested by Magnus [0] 


S = e 




( 11 ) 


Here the Hermitian operator F{t) can be represented as a series of multiple commuta¬ 
tors with t-integrations 


m = v{t) - Y'(t),v{t)] - li v(t).|r(t) ,y(i)]| - h ir(i).r(0] .r(t)| 

+4| r(t).r(t),r(i)] ,v'(t)]] 

+ 4l |r(t), i/(i)]| , i/(i)] + 4[ ir(i), v{t)], y(t)|| +... ( 12 ) 

One important advantage of this representation is that the A-operator (ED is manifestly 
unitary. It follows from equations © and ED that operators E(f) and F{t) are related 
to each other 
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m 




( 13 ) 


so finding F{t) or S(t) are eqnivalent tasks. 

The S'-operator and the Hamiltonian provide two different ways to describe dy¬ 
namics. The S'-operator represents only “integrated” time evolntion from the remote 
past to the distant fntnre. The knowledge of the S'-operator is snfficient to calcnlate 
the scattering cross-sections as well as energies and lifetimes of stable and metastable 
bonnd states.^ However, in order to describe the time evolntion and the wavefnnctions 
of bound states the full interacting Hamiltonian H is required. 

It can be shown [7] that two Hamiltonians H and H' related to each other by a 
unitary transformation 


H' = 

yield the same scattering S" = S' as long as condition 


lim = 0 (14) 

t^±oc ^ ’ 

is satished. Such Hamiltonians H and H' are called scattering-equivalent. The energy 
spectra of two scattering equivalent Hamiltonians are identical. However, the eigenvec¬ 
tors are different and, according to eq. dH), the corresponding descriptions of dynamics 
are different as well. Therefore scattering-equivalent theories may be not physically 
equivalent. 

Calculations of bound states, time evolution and scattering is a routine practice in 
non-relativistic quantum mechanics. However, the situation is less certain in relativistic 
quantum held theories. As mentioned in Introduction, the Hamiltonian of renormalized 
QFT is inhnite. Therefore, it is not immediately clear if the above formulas O - (HI 
remain valid for the high energy relativistic phenomena, and what modihcations, if 
any, should be introduced in quantum theory to take into account the variable number 
of particles. In this paper we illustrate the difficulties encountered in quantum held 
theories by analyzing a simple model theory with variable number of particles. We 
will demonstrate that the Hamiltonian can be redehned so that there is no need for 
renormalization and usual quantum mechanical techniques remain applicable even in 
the relativistic case. 

^The latter two quantities are represented by positions of poles of the S'-operator on the complex 
energy plane. 
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3 Model theory 

Our model theory describes two kinds of particles. These are massive spinless fermions 
which will be called electrons and massless bosons with zero helicity, which will be 
called photons. Here we disregard the spin and polarization degrees of freedom as they 
are not so important for the discussion of renormalization. To allow for creation and 
annihilation of particles, the system is described in the Fock space which is built as 
a direct sum of sectors with various numbers of particles. For example, if we denote 
|0) the no-particle vacuum state, Hei the one-electron Hilbert space and Hph the one- 
photon Hilbert space, then the Fock space can be written as an inhnite direct sum 


n = \0)®nel®nph®(nel®nph)(B(nel(^asymnel)(B(nph(B)symnph)...{l5) 

The anticommutation and commutation relations for particle creation and anni¬ 
hilation operators (a^, Up for electrons and cj,, Cp for photons, respectively) are, as 
usual. 


{Op, Up/} 

= <^(p-p') 

(16) 

[ck,4,] 

= <^(p-p') 

(17) 


= {ai’®P'} = 0 

(18) 

[Cki Ck^] 

[cL 4'] = 0 

(19) 


[aj>) 4] = [ap> Ck] = [ftp, 4] = [ttp, Ck] = 0. 

The full Hamiltonian H = Hq + Vi is the sum of the free Hamiltonian 


Ho = J dpupolap + j dk|k|cj,Ck (20) 

kT^O 

(where Up = ^ind |k| are one-particle energies of electrons and photons, 

respectively) and the interaction, which we choose in the following form 


Hi = e(27r) f ^^^aj,4ap+k + e(27r) 3/2 f 

V|k| 4o 


kT^O 


dpdh 

Tw 


®p®p—kCk 


( 21 ) 
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The coupling constant e is the absolute value of the electron charge. Here and in what 
follows the perturbation order of an operator (= the power of the coupling constant 
e) is shown by the subscript. For example, the free Hamiltonian Hq does not depend 
on e, so it is of zero perturbation order; Vi is of the first perturbation order, etc. The 
number of electrons is conserved by the interaction m, but the number of photons 
is not conserved. So, this theory is capable of describing important processes of the 
emission and absorption of photons. 

In this paper we will discuss ultraviolet divergences associated with the interaction 
m- However, we will skip completely the discussion of “infrared divergences” which 
are related to the zero mass of photons and singularities in m- The easiest 

way to avoid infrared problems in practical calculations is to assign a small non-zero 
mass to photons. 

In section El we are going to calculate the scattering operator Q with the above 
Hamiltonian. Before doing that, some remarks are in order. The operator S is obtained 
as a sum of products (HDD (or commutators (HI) of interactions Vi(t) = e 
with f-integrations. Each term in these expressions can be written as a normally 
ordered product of N creation operators and M annihilation operators a.^ 

VNMit) = J [dq]T)ArM[q]e**^'^'^(‘^(h(P7VM[q])aq/^ ■ • • • • • «qM (22) 

where integration is carried over momenta of all created and annihilated particles [dq]. 
The momentum conservation law is guaranteed by the delta function in ( 1221 ) whose 
argument is the sum of momenta of created particles minus the sum of momenta of 
annihilated particles 


N M 

PiVM(q'i,..., q)v, qi, • • •, qv/) = Z)~ H (23) 

i=l j=l 

Usually, we will perform explicit integration over one momentum q'^ which removes the 
delta function and expresses q'^ as a linear function of other momenta. The argument 
of the exponent in (122|l contains the energy function 

N M 

^iVM(q'i,..., q)v, qi, • • •, qv/) = Z ^q' “ Z ^q. (24) 

i=l j=l 

which is the difference of energies of particles created and annihilated by U/vm- 
is a numerical coefficient function. 


^Here symbols and a refer to generic creation and annihilation operators without specifying the 
type of the particle. The pair of integers {N,M) will be referred to as the index of the term Vnm- 



Suppose that a term V)vm has coefficient function D^m, then we introduce a useful 
notation V)vm°C for fh® operator whose coefficient function is a product of D^m 
and a function ( of the same arguments 


D'nmH = -DiVM[q]C[q] 

Then, a t-dependent regular term V^Mif) can be written as 


VNuit) 


and its dehnite t-integral is 




Vnmo 


VNM{t) 



VNM{t)dt 


2tiVnm o S{Enm) 


(25) 


Eq. (j25jl means that each term in V^Mit) is non-zero only on the hypersurface of 
solutions of the equation 


Enm{(Ii, • • •) qw) qi) • • •) qAf) — o 

if such solutions exist. This hypersurface in the momentum space is called the energy 
shell of the term VNMif)- Note that the scattering operator (jHl) is non-trivial only on 
the energy shell, i.e., where the energy conservation condition holds. 


4 Three types of operators in the Fock space. 

In this section we would like to get a further insight into the nature of operators in the 
Fock space by dividing them into three groups depending on their index {N,M). We 
will call these types of operators renorm, phys, and unphys. 

Renorm operators have either index (0,0) (a numerical constant C) or index 
(1,1) in which case the same type of particle is created and annihilated. The most 
general form of a renorm operator in our theory is^ 

^Here we write just the operator structure of R omitting all numerical factors, indices, integration 
and summation signs. Note also that the interaction operator eu commutes with the operator 
of charge Q = —e f dpa^Up, so any product or commutator of terms derived from Vi should also 
commute with Q. This is not true for terms like a^c and c^a, so they are not allowed in the theory. 
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R 


a^a + c'*'c + (7 


(26) 


The free Hamiltonian dm) is an example of a renorm operator. The class of renorm 
operators is characterized by the property that the energy function is identically 
zero. So, renorm operators always have an non-empty energy shell where they do not 
vanish. Regular renorm operators do not depend on t. 

Phys operators have at least two creation operators and at least two destruction 
operators (index {N,M) with N > 2 and M > 2). In this case the energy shell is 
non-empty. For example, the energy shell for phys operator Op+k^q-k^p'^q is given by 
the set of solutions of equation Wp+k + <^q-k = + i^q which is not empty. 

Unphys operators have index (1,A^ > 2) or (A^ > 2,1). It can be shown that 
the energy shell is empty for unphys operators. For example, the interactions m are 
unphys. The energy shell equation for these two terms ujp + |k| = (x^p+k does not have 
a solution, because there are no photon states with zero momentum. 

Renorm, phys, and unphys operators exhaust all possibilities in our theory, there¬ 
fore any regular operator V must have a unique decomposition 


V{t) = y’" + I/“"^’(f) + VP^{t) 

The rules for calculations of commutators, derivatives and f-integrals with different 
operator types are summarized in Table 9.2. For example, the f-integrals of phys and 
unphys operators are given by formula 


V{t) = V{t)o^ 


(27) 


Table 1: Operations with regular operators in 
U=unphys, R=renorm, NR=non-regular.) 

the Fock 

space. 

(Notation: 

P=phys, 

Type of operator 

A 

\A,P] 

\A,U] 


dt 

A 


P 

P 

P+U 

P 

P 

P 

P 

U 

P+U 

P+U+R 

U 

u 

U 

0 

R 

p 

U 

R 

0 

NR 

oo 
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+1 


(a) 


+ 1 


(b) 



Figure 1: Diagram representation of two terms in the interaction operator Vi (eq. 

(EU). 

5 Diagrams in the model theory. 

Our goal in this section is to introduce the diagram technique which greatly facilitates 
perturbative calculations of scattering operators m and m- Let us graphically rep¬ 
resent each term in the interaction operator m as a vertex (see Fig. 1). Each particle 
operator in m is represented as an oriented line or arrow. The line corresponding to 
the annihilation operator enters the vertex, and the line corresponding to the creation 
operator leaves the vertex. Electron lines are shown by full arcs and photon lines are 
shown by broken arrows. Each line is marked by the momentum label of the corre¬ 
sponding particle operator. Free ends of the electron lines are attached to the vertical 
“order bar” on the left hand side of the diagram. The order of these external lines 
(from bottom to top of the diagram) corresponds to the order of particle operators in 
the interaction term (from right to left). An additional numerical factor is indicated 
in the upper left corner of the diagram. 

The f-integral Vi{t) differs from Vi{t) only by the factor iEy^ (see eq. (|2Z|)) which 
is represented in the diagram by drawing a box that crosses all external lines. A line 
entering (leaving) the box contributes its energy with the negative (positive) sign to 
the energy function Ey^- The diagram representation of the integral 


vAt) 


is shown in Fig. 2. 


ie 

(27r)3/2 

ie 

(27r)3/2 


dpdk . t 

dpdk e-*h^P-|k|-<^p-k) 


p-k 




( 28 ) 
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+1 


(a) 


+ 1 


(b) 



Figure 2: t-integral Vi{t) 


The product of two operators AB is represented by simply placing diagram B below 
diagram A and attaching the external electron lines of both diagrams to the same order 
bar. For example, the diagram for the product of the second term in ()21j) (Fig. 1(b)) 
and the hrst term in (OK|l (Fig. 2(a)) 


ViVi oc (aJ,ap_kCk)(o]j4,aq+k') + • • • (29) 

is shown in Fig. 3(a). This product should be further converted to the normal form, i.e., 
all incoming lines should be positioned below the outgoing lines. Due to the relations 
(uni) - (uni) each exchange of positions of the electron particle operators (full external 
lines on the diagram) changes the total sign of the expression. Each permutation of 
annihilation and creation operators (incoming and outgoing lines) of similar particles 
creates an additional expression and a new diagram in which the swapped lines are 
joined together. Using these rules we hrst move the photon operators in to the 
rightmost positions, move the operator to the leftmost position, and add another 
term due to the anticommutator {ap_k, = 5(q — p + k). 


Ui\h oc a|ja{,ap_kaq+k'Ck4, + (j(q - p + k)a{,aq+k'Ck4/ 

ClqQ.pQ.p_kOq+k'0kCk' T ®p®p—k+k'^kCk' T . . . (30) 

This expression is represented by two diagrams 3(b) and 3(c). In the diagram 3(b) the 
electron line marked q has been moved to the top of the order bar. In the diagram 
3(c) the product h(q — p + k) and the integration by q are represented by merging 
or pairing the incoming electron line carrying momentum p — k with the outgoing 
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Figure 3: Normal product of operators in Fig. 1(b) and 2(a). 


electron line carrying momentum q. This produces the internal electron line carrying 
momentum p — k between two vertices. 

In the expression jsni), the electron operators are in the normal order, however, 
the photon operators are not. The next step is to bring the photon operators to the 
normal order 


kj lh 0^ ®qf^p®p~k®q+k'Cj^'Ck T ®q®p®p~k^^q+k'^(k k) 

+ Opap~k+k'4'Ck + oJ,ap~-k+k'5(k' - k) + ... 

OqClpCip—k®q+k'Cj^'Ck T ®q®p®p—k^^q+k 
T Q-p^p—k+k'Ck'Ck T Q-p^p T . . . 

According to equation m, the normal ordering of photon operators in 3(b) yields 
diagrams 3(d) and 3(e).® Diagrams 3(f) and 3(g) are obtained from 3(c) in a similar 
manner. 

Using diagrams, with some practice, one can perform calculations of scattering 
operators m and m much easier than in the usual algebraic way. During these 

®We relabel the external momenta in Fig. 3(e) for future convenience. 
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diagram manipulations we, actually, do not need to keep momentum labels of lines. 

The algebraic expression of the result can be easily restored from an unlabeled diagram 

by following these steps: 

(I) Assign a distinct momentum label to each external line, except one, whose mo¬ 
mentum is obtained from the momentum conservation condition. 

(II) Assign momentum labels to internal lines so that the momentum conservation 
law is satisfied at each vertex. If there are loops, one needs to introduce additional 
independent loop momenta^ 

(III) Read external lines from top to bottom of the order bar and write corresponding 
particle operators from left to right. 

(IV) For each box, write a factor i{Ef — Ei)~^, where Ef is the sum of energies of 
particles going out of the box and Ei is the sum of energies of particles coming 
into the box. 

(V) Write a factor where Ey is the energy function of the diagram which is 

the sum of energies of all outgoing external lines minus the sum of energies of all 
incoming external lines. 

(VI) For each vertex introduce a factor , ^ -, where k is the momentum of the 

photon line attached to this vertex. 

(VII) Integrate the obtained expression by all independent external and loop momenta. 

6 Electron-electron scattering. 

Let us now try to extract some physical information from the above theory. We will 

calculate low order terms in the perturbation expansion m for the F'-operator. 


Si(t) = zVi(t) (31) 

^2it) = -iVi{t)Vi{t)rP - (32) 

To obtain the corresponding contributions to the F-operator we need to take t-integrals 
of these expressions 




1 -1- Si(t) -1- S2(f) -|- ... 


^see diagram 3(g) in which k is the loop momentum. 
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Note that the right hand side of m and the hrst term on the right hand side of 
(El are nnphys. Their energy shell is empty, so, according to Table 1, they do not 
contribute to the S'-operator. 

Operator on the right hand side of eq. dSl has two terms corresponding 

to two types of scattering processes allowed in the 2nd perturbation order.® The term 
of the type a'^c^ac (see, e.g., fig. 3(f)) annihilates an electron and a photon in the 
initial state and recreates them (with different momenta) in the final state. So, this 
term describes the electron-photon (Compton) scattering. In this paper we will focus 
on the other term which describes the electron-electron scattering. Let us consider in 
more detail the second-order contribution to this process (see fig. 3(e)) 


S'aia^a^aa] = 27 r y dp(iqdkh(a;p_k + ^^q+k - - a;p)T>2(p, q, k)a^„j^a^^j,«q«P 

The coefficient function in dSl can be read from the diagram, according to the rules 

(I) - (VII), 


A(p,q,k) 


ie^ 1 

(27r)3 |k|(|k| -hcup-k - lOp) 


(33) 


In the non-relativistic approximation (p, q,k me), the coefficient function D 2 {p, q, k) 
has singularity |k|~^ which is characteristic for scattering of two electrons interacting 
via repulsive Coulomb potential 


47r|ri - r2| 

So, our model theory is quite realistic. 

7 Mass renormalization. 

Next consider the third term on the right hand side of eq. ()32|1 . It is given by the 
diagram in fig. 3(g). According to rules (I) - (VII) this diagram is represented by the 
expression 




le 


( 21 ,) 


(ipdk 




(^p—k ^p “1“ 


(34) 


®We do not need to consider unconnected diagrams, like fig. 3(d), because they describe two or 
more disjoint scattering processes. 
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There are serious problems with this term. First, the loop integral by k is divergent 
because the integrand in (P|) has asymptotic behavior oc /c ^ at large k. However, even 
if the integral were convergent, the presence of a renorm contribution in the operator 
S(f) is unacceptable, because, according to Table 1, the f-integral of any renorm term 
is infinite. So, the scattering phase in the second order E 2 is infinite, which is absurd. 

Moreover, if we continued calculations dSH) - (IS21) to higher perturbation orders we 
would find out that even phys terms in S(f) become inhnite due to divergent loop 
integrals, so the theory with the Hamiltonian (PU|) - (EH) is seriously flawed. 

In a consistent theory we must require that 


E""’" = 0 (35) 

Therefore, operator must be purely phys 

^ ^ (36) 

It was shown in ref. PP that condition EBD is equivalent to requiring that the F-operator 
leaves the vacuum and one-particle states invariant 


^| 0 ) = 

|0) 

(37) 

Sal\0) = 

oJ.IO) 

(38) 

^ 4 | 0 ) = 

4lo> 

(39) 


This is the mass renormalization condition of the traditional renormalization theory. 
To satisfy this condition, we must modify the Hamiltonian EHD- EH) by adding certain 
unphys U and renorm R counterterms to the interaction operator Vi.® In other words, 
we are saying that the original Hamiltonian 


H = H^ + Vi 


(40) 


is not correct, and the modihed Hamiltonian with renormalization counterterms 


= Ho + Vi + U + R (41) 

®It appears that in our model theory the renormalization is achieved by unphys and renorm coun¬ 
terterms only. In the general case, e.g., in QED, phys counterterms should be added as well. 
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better describes interactions between particles. To comply with eq. dsni) , we must 
choose the counterterms in such a way that operator 


E"(t) = iV\t) - V\t) V\t) + ■ ■ ■ (42) 

does not contain renorm terms. From Table 1, it is clear that renorm terms in 
may appear due to the presence of renorm and unphys terms in V‘^ and their products. 
So, in order to satisfy eq. there should be such a balance between unphysical 

and renorm terms in V^lt) that all renorm terms in cancel out in all orders 

of the perturbation theory. The mass renormalization is achieved by adding renorm 
counterterms in even orders: R 2 , R 4 , etc. The charge renormalization procedure will 
be discussed in section |H1 It requires addition of unphys counterterms in odd orders 
Usit), U^it), etc. We take these considerations into account by writing the general 
expression for the Hamiltonian of the renormalized theory^° 


H%t) = Ho + V^it) 

where 


V%t) = V,{t) + R 2 + U3{t) + R4 + ... (43) 

We obtain formulas for S^(f) {i = 1, 2, 3,...) by inserting interaction in (ITTH) and 
collecting terms of equal order. 


E5(f) 

= ^V7(^) 

(44) 


= + tR2 

(45) 


= -tW(t)W(t)W(f) - i?2W(t) - Hl(t)^ + z[/3(t) 

(46) 

m) 

= a4{t)-U3{m{t)-V4{t)Us{t) + iR4 

(47) 

where we denoted 




= Hi(f)l/i(t) W(t)W(f) -^W(t) l/i(f)^ 

- tVi{t) R2Vi{t) -tR2 Vi{t)Vi{t) , 

course, we are looking for a f-independent Hamiltonian H°. Although, at intermediate calcu¬ 
lation steps it is convenient to keep all operators t-dependent, as in eq. ©, in the end we should set 
t = 0. 
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Now we go order-by-order and choose counterterms • • • so that renorm terms 

are eliminated from the left hand sides of eqs dSD - (IIZD- The first-order term is 
unphys, so there is no need for renormalization in the first order. To ensure that 
does not have a renorm part we choose the counterterm 


i?2 = 

(see diagram 3(g)). With this choice, we can rewrite the contributions to the iS-operator 
in the 4 lowest orders 

gw = 0 

gw = 

gw = -i vmvimm”*" ( 48 ) 

ES(() = ai(t)-U3(t)Vy)-\\{t)U3H)+iRi, (49) 


where 


= ^i(f)V)(t)(W(t)W(t))^+" + W(t)W(t)(W(t)W(t))^^" 

V y V y 

V V 

- iV^{t)R2V^(t) -iR2V^(t)V^(t) -iV^(t)Vi{t)^ 

v ~ v ~ ^ ^ ^ 

= W(t)W(t)(W(t)W(f))^+"-^W(t)/^2W(t) 

'- V- '---' 

- iR2Vi{t)Vi{t) (50) 


and the superscript p R u denotes the sum of phys and unphys terms. The term on 
the right hand side of (1481) has odd number of particle operators, hence it is free of 
renorm parts (which have either zero or two particle operators). Therefore, no renorm 
counterterms should be added there. Just as we did in the 2nd order, we can choose a 
renorm counterterm i ?4 which simply cancels all renorm terms which may be present 
in the first three terms on the right hand side of ()49j) . With the above choices, the 
operator S‘^(f) does not contain renorm terms up to the 4th order, as required by the 
mass renormalization condition,and the expression for E 4 simplifies 

^^Our analysis demonstrates that that mass renormalization is always necessary when interaction 


18 








































' -V-' 

mt){vYt)vYt)r^vYt)Y'^-{UsmiiYr’^-ivYm^^ ( 51 ) 

'-V-' '-V-" '-V-" 

Using the diagram technique we find that the a^a^aa part of the first two terms on 
the right hand side of eq. (ED is represented by 7 diagrams shown in Fig. 4. From 
diagram rules (I) - (VII), we obtain the contribution of diagrams 4(a) - 4(b) to the 
coefficient function D for the electron-electron scattering on the energy shell. 



D4 '’^(p,q,k) = 


^e 


(27r)6(a;p_k - + k)k 


dh 1 


+ 


EF' 


(62) 


where B = a;p_h — oJp + h, C = Wp^k — <^p-h-k — h, E = Wq-h — a;q -|- h, and 
F = Wq+k — <^q+k-h “ h. Unfortunately, this contribution is infinite (at large values 
of h the integrand behaves as h~Y- Thus we conclude that the mass renormalization 
procedure described above has not removed all divergences. To solve this problem we 
need to perform the second renormalization step known as the charge renormalization 
procedure. This step is explained in the next section. 


8 Charge renormalization 

First note that the divergent terms (HI have a singularity k ^ at fc —0. As we know 
from section El such a singularity is responsible for the low-energy electron-electron 
scattering at large distances. From classical physics we also know that long-distance 
interactions between charged particles depend on (in our language, they are of the 
second perturbation order) and they are accurately described by the 2nd order term 
dSSD- Non-zero terms D^Y mean that the charge of the electron is modified by the 
interaction. Actually, this modification is infinite, but even finite values of the terms 
like are inconsistent with the classical limit. So, we will postulate that in orders 
higher that 2nd, singular coefficient functions like (jn2D should not be present at all, 
whether they are infinite or finite. This is the charge renormalization eondition. To 
eliminate the divergent contribution ()52ll we can choose unphys 3rd order counterterms 

contains unphys terms, like trilinear operators in eq. ED and in more realistic theories, such as QED, 
QCD, and Standard Model. Products or commutators of such terms in m or m give rise to renorm 
terms in the S'-operator which should be compensated by adding renorm terms to the Hamiltonian. 
The operator structure of these terms is the same as in Hq (EOI). This results in the difference of 
masses of bare and dressed particles, as discussed at the end of section [S] 


19 


















1 


1 



Figure 4: Contributions with the operator structure a^a^aa to the hrst two terms on 
the right hand side of eq. m- 
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Figure 5: Charge renormalization counterterms Uslt). 


U 3 . Two diagrams describing these counterterms are shown in Fig. 5(a) and 5(b). As 
expected, they are inhnite, and expressions 


in (im exactly cancel unwanted inhnite diagrams 4(a) and 4(b) on the energy shell. 

Surviving diagrams 4(c) - 4(g) are the 4th order radiative corrections to the electron- 
electron scattering. On the energy shell they yield the following contribution to the 
coefficient function of S(f) 


D'f ®^(p,q,k) 
where 


A = 

- <^q-h 

— 

^^q + 

h 




B = 

- ^p~k 

-h 

Co^p 

+ 

|h 

+ 

k| 

C = 

= l^q-h 

+ 

k 

+ 

|h 

+ 

1 

1 

D = 

“ ^^p+h 

+ 


— 

CUq 


UJ-p 

E = 

= l^q+k 

— 

<^q + 

|h 

+ 

k| 

+ h 

F = 

= l^q+k 

+ 

l^p-h- 

-k 

+ 

h- 

— Co^q Up 

G = 

= ^q+k 

+ 

^^p+h 

+ 

|h 

+ 

S — CUq — 


ie^ r dh 1 1 1 1 1 1 

(27r)6 J h\h + k\A^^ ^ ^ ^ ^ ^ ^ ^ 
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The integrand in dSSl) is proportional to h and the integral is convergent at large 
values of the loop momentum h. Therefore is hnite on the energy shell. So, 

in the renormalized theory the electron-electron scattering is represented by a finite 
iS-operator which up to the 4th order is described by the coefficient function 


D — Z ?2 + . 

where D 2 is given by eq. O and D 4 is given by eq. (ESD. 

The renormalization technique presented above has not been applied to realistic 
theories, such as QED, yet. However, the close analogy between renormalization issues 
encountered in our toy model and in QED allows us to speculate that similar renor¬ 
malization steps can be repeated in all perturbation orders in QED, so that hnite and 
accurate values for various scattering amplitudes can be obtained. Note that in our 
approach the counterterms in the Hamiltonian BD have the same operator structure 
as the terms in the original Hamiltonian (HOD. In relativistic renormalizable theories, 
like QED, even stronger statements can be made: the mass renormalization countert¬ 
erms differ from the one-particle energy terms in Hq only by a constant (but inhnite) 
factor,and the charge renormalization counterterms are equal to the interaction term 
V in the original Hamiltonian multiplied by a constant inhnite factor. This conclusion 
is true in all perturbation orders. Thus, if the original Hamiltonian depends on hnite 
masses mi, m 2 ,... and coupling constants (charges) ei, 62 ,... 


H(mi,m2,... ;ei,e2,...) (54) 

Then the addition of renormalization counterterms is equivalent to simply modifying 
(making them inhnite) the values of masses and charges. The Hamiltonian after 
renormalization has the same functional form as H where parameters mi, m 2 ,... and 
Cl, 62 ,... substituted by renormalized (inhnite) values rhi, m 2 ,... and Ci, 62 ,... 


if'" = i/(mi, m2,...; 61,62,...) ( 55 ) 

The traditional approach ohers the following physical interpretation of these results. 
The creation and annihilation operators present in the theory describe so-called bare 
particles having inhnite masses and charges. The interaction between bare particles 
drastically change their properties. For example, bare electrons constantly emit and re¬ 
absorb virtual photons and electron-positron pairs. So, bare electrons are surrounded 
by a coat of virtual particles. The mass of this coat is inhnite, thus compensating 
the inhnite mass of the bare particle. The virtual particles in the coat shield the true 

^^From this point of view the mass renormalization was discussed in ref. jH). 
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(infinite) charge of the bare electron. The resulting mass and charge of the dressed 
particle is hnite and exactly equal to the parameters measured in experiment. In 
experiment, we never see the bare particles and their virtual coats, we only see the 
dressed particles. 

The Hamiltonian is formally inhnite, but these inhnities cancel when the S- 
operator is calculated by formula 0, and accurate results are obtained for scattering 
amplitudes and energies of bound states. However, the inhnities do not cancel when 
one tries to hnd the eigenvectors of the Hamiltonian or calculate the time evolution 
operator exp{iHH). Direct application of this formula would lead to unphysical results. 
For example, if we calculated the time evolution of the simplest one-electron state we 
would obtain that this states dissociates into a complex linear combination of states 
over time 


e*^V|0> = {l + iHH+...)a^\0> 

oc (1 -f it{a^a + c^c + a^c^a + a^ac) -|- .. .)a^|0 > 

oc a^|0 >-|-a'*'c'*'|0 >-|-a'*'c'^c^|0 > -I- ... (56) 

Moreover the coefficients multiplying the terms in (jini) are given by divergent integrals. 
Therefore, with the inhnite Hamiltonian the basic formulas of quantum theory 
o and o become useless for practical calculations. However, in most cases this 
is not a cause of trouble. As discussed in section [2l most experiments in high energy 
physics are concerned either with bound state energies or with scattering cross-sections. 
To calculate these properties, the knowledge of the S'-operator is sufficient and the 
renormalized theory works hne. For approximate description of the time evolution in 
low energy, e.g., atomic, systems, one can use heuristic approaches, such as those based 
on the Dirac-Fock Hamiltonian. 


9 Construction of the dressed particle Hamiltonian 

We see two fundamental problems with the traditional renormalization theory pre¬ 
sented above. First, from the practical point of view, the absence of a well-dehned 
hnite Hamiltonian does not allow one to study the time evolution of interacting states. 
One can expect that with advancement of experimental tools, the time-dependent in¬ 
formation from the region of interaction will soon become available. However, accurate 
analysis of this information is not possible with ill-dehned Hamiltonians of renormalized 
quantum held theories. Second, from the theoretical point of view, the traditional 

^^The way to remove the ultraviolet divergences from interaction operators, or to “renormalize” 
the Hamiltonian, was suggested by Glazek and Wilson in their similarity renormalization approach 
P]. However, their Hamiltonians still have unphys terms |H], and result in “instability” of vacuum 
and one-particle states similar to that shown in eq. 
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approach operates with notions of bare and virtual particles that are, in principle, 
non-observable. It is desirable to have a theory formulated directly in terms of dressed 
particles and their interactions without mentioning the non-observable bare and vir¬ 
tual particles at all. This idea was hrst realized in the dressed particle formalism by 
Greenberg and Schweber p. 

In the rest of this paper we will discuss the RQD approach PUP which combines 
the similarity renormalization and dressed particle ideas. RQD is capable to hx the 
both problems of the traditional renormalization formalism mentioned above. Let us 
hrst focus on the dressing part of RQD. There are two approaches to dressing with 
different interpretations but equivalent physical results. One approach ma tries to 
hnd an explicit unitary (dressing) transformation connecting creation and annihilation 
operators of bare particles (a^, a, c^, and c) with creation and annihilation operators of 
dressed particles C*!, and C) and hnd a function / which expresses the Hamil¬ 

tonian of the renormalized theory through dressed operators 


H^ = f(A\A,C\C) (57) 

Another approach PPP is to pretend that operators a\ a, c\ and c already describe 
the dressed particles and try to hnd a new hnite Hamiltonian whose functional 
dependence on bare particle operators is given by the same function / which expresses 
the dependence of on the dressed operators 


= f{a\a,c\c) (58) 

We will stick to the second approach in this paper. There are two ways to proceed. 
First, we can try to hnd H'^ by applying a unitary dressing transformation to H^. This 
way was described in ref. [3]. In this paper we are choosing another route: we will 
simply ht the Hamiltonian H'^ to the hnite scattering operator S'^ known from the 
renormalized theory. 

There are three requirements that we want to satisfy when looking for the “dressed 
particle” Hamiltonian = Hq + = Hq + V 2 + + Vf + .... 

(A) is scattering-equivalent to H'^. 

(B) is finite. 

(C) is phys. 

The condition (A) is understandable as we know that the R-operator of renormal¬ 
ized theories (QED, Standard Model) agrees well with experiments, and we would like 
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to preserve this agreement in a theory with the new Hamiltonian H’^. Using eqs. m 
and m we can write 


-ilogS' = £JW+£3W+£4W + ^-- 

-iiogS" = gw Vj'M v-y;)] ... 

Then we see that condition (A) implies the following infinite set of relations between 
Vf{t) and F^{t) on the energy shelh'^ 


vf{t) 


Fnt) + Q^{t), z>4 


(59) 

(60) 

(61) 

(62) 


where Qi(t) denotes a sum of multiple commutators of from lower orders (2 < 

j < i — 2) with t-integrations. We have expressions for operators ^^(t) and ^^(f) in the 
renormalized theory (eqs. (jddj) and (jSSl))- Clearly, we can perform similar calculations 
for operators F^{t) whose values on the energy shell are present on the right hand sides 
of eqs. (Ei - Alternatively, we can express F^{t) through E'^(t) using relationship 

m- 

Operators F^{t) are, of course, finite. This immediately implies that V^it) and 
V^{t) are finite on the energy shell. Moreover, from eqs. and it follows 

that F 2 {t) and F^{t) are phys on the energy shell. Therefore, V^it) and V^{t) can 

be also chosen phys on the energy shell. This proves that conditions (B) and (C) are 
satisfied on the energy shell up to the 3rd order. However, eqs. (Pi) - (in2i) tell us 
nothing about the behavior of V^it) and V^{t) off the energy shell.The same is 
true for interactions in higher orders: the interaction operators V^{t) off the energy 
shell remain undetermined by our condition (A). This freedom in choosing interactions 
simply means that different choices of the coefficient functions of Vf'{t) off the energy 

use of the Magnus expansion for deriving V/{t) is preferable to the old-fashioned or Dyson’s 
formulas, because the right hand sides of eqs. PI - 113) are expressed through commutators, so 
they are manifestly Hermitian. In addition, this approach automatically generates cluster separable 
interactions Vf{t) (see ref. |1]). 

^^Note that in applications we are dealing primarily with interactions near the energy shell, because 
in most processes, excluding very short virtual events, the total energy stays unchanged at all times. 
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shell result in scattering equivalent Hamiltonians. Since our experimental knowledge 
about interactions in high energy physics is limited to their effects on scattering, we are 
allowed to freely choose the behavior of Vf{t) off the energy shell without any danger 
to get into contradiction with experiment. This freedom is exactly what is needed to 
satisfy conditions (B) and (C) in all orders. As we will see shortly, it is important 
to choose the behavior of Vf'{t) such that their coefficient functions fall off rapidly 
when the arguments move away from the energy shell.For example, we can choose 
the coefficient functions being proportional to Q = exp(—where 7 is a positive 
constant and E is the energy function. Then the electron-electron interaction in the 
second order takes the form 


= -;KWViW'‘°C2 


(2ir) 


dpdqdk- 


^-■yE^ ^-itE 


k|(|k| -Fcup-k - ojp) 


®p-k®q+k®q®P 


(63) 


where E = Wp^k + <^q+k — i^q — is the energy function. 

There are no terms of the type a^a^aa in the 3rd order term F^{t). Therefore the 
lowest (4th) order radiative correction to the electron-electron interaction should 
be obtained from eq. (inn). As discussed above, the value of T|(t) on the energy shell 
is hnite. Let us now consider the term 


(64) 


on the right hand side of eq. (jnn). First, we note that, according to Table 1, this 
commutator is phys. Second, as we agreed above, the coefficient functions of 
fall off rapidly outside the energy shell. Without this condition, the loop integrals 
encountered in calculations of (EH) could be divergent. However, in our case these 
integrals are convergent: when the loop integration momentum goes to infinity, the 
V-fit) factors in the integrand go away from the energy shell, i.e., rapidly fall off. 
This guarantees that interaction V/(t) is phys and finite on the energy shell. Its 
coefficient function off the energy shell should be again chosen to decay rapidly to 
ensure the convergence of loop integrals in higher order operators Qi{t), where V/(t) 

the representation where Hq is diagonal, this condition is equivalent to bringing the matrix of 
the Hamiltonian to the band-diagonal form, which is the central idea of the similarity renormalization 
method P). 

^^Note that we should be careful not to set 7 to infinity. In this case, the coefficient functions of 
interaction operators Vf' become non-differentiable. This means that interaction is not separable El, 
and the scattering theory formalism from section PI is no longer applicable. 
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may contribute. These arguments can be repeated in higher orders, which proves that 
the dressed particle Hamiltonian H‘^ is free of ultraviolet divergences. 

The contribution on the right hand side of eq. iH) is well-defined near 

the energy shell, but this is not true for the contribution This commutator 

depends on the behavior of V 2 everywhere in the momentum space. So, it substantially 
depends on our choice of ^2 off the energy shell. There is a great freedom in this choice 
which is reflected in the uncertainty of even on the energy shell. It is interesting to 
note that although the off-shell behavior of the 2nd order interaction and the on-shell 
behavior of the 4th order interaction cannot be separately determined in our theory, 
they are connected to each other in such a way that the ambiguity of the interaction 
does not affect the S'-matrix for the electron-electron scattering. 

The above construction does not allow us to obtain full information about V^-. 
The off-shell behavior of interactions is rather arbitrary, and the on-shell behavior 
can be determined only for lowest order terms. However, this uncertainty is perfectly 
understandable: It simply reflects the one-to-many correspondence between the S- 
operator and Hamiltonians (see section ED. It means that there is a class of finite 
phys interactions all of which satisfy our requirements (A) - (C) and can be used 
for S'-matrix calculations without encountering divergent integrals. In order to find 
the unique Hamiltonian correctly describing the dynamics of particles, the theoretical 
predictions should be compared with time-resolved experimental data. However, our 
model theory is not sufficiently accurate to be comparable with experiment, and further 
efforts in this direction require building the dressed particle formulation of the full¬ 
blown quantum electrodynamics 13 El- 


Table 2: Examples of interaction terms in the dressed particle Hamiltonian dSHD- Bold 
numbers in the third column indicate perturbation orders in which the operators can 
be unambiguously obtained near the energy shell as discussed in section El 


Operator 

Physical meaning 

Perturbation 

Orders 

a^a^aa 

Elastic potentials 

electron-electron 

2,4,6,... 

aY^ac 

electron-photon (Compton) 

2,4,6,... 

a^a^a^aaa 

3-electron potential 

4,6,... 

a^aY^aa 

Inelastic potentials 

bremsstrahlung in electron-electron collisions 

3,5,... 

a^a^aac 

photon absorption in electron-electron collisions 

3,5,... 


In contrast to the original Hamiltonian (HOD, there seems to be no way to write 
the dressed particle Hamiltonian H‘^ in a closed form. From the derivation outlined 
above it is clear that in higher perturbation orders there are more and more terms with 


27 





increasing complexity in the interaction operator 


= Ho + V^ 

= Hq + a^a^aa + a)c^ac + o)o)a)aaa + o)a)aac + a)o)c^aa + ... (65) 

Some of them are shown in Table 2. However, all these high order terms directly reflect 
real interactions and processes observable in nature. For example, the term a)a^c^aa 
{bremsstrahlung) describes creation of a photon in electron-electron collisions. In the 
language of classical electrodynamics, this can be interpreted as radiation due to the 
acceleration during interaction of charged particles and is often referred to as the 
radiation reaction force. The Hermitian-conjugated term a^a^aac describes absorption 
of a photon by a colliding pair of charged particles. 

Note that in contrast to the traditional renormalization approach in which the mass 
renormalization condition (jai was ensured by maintaining a balance between unphys 
and renorm terms in the interaction, in the RQD Hamiltonian these terms are absent 
altogether. The interaction operator is purely phys which guarantees that scattering 
operators and F'^ are phys too. The phys character of the interaction means that 
it yields zero when acting on the vacuum and one-particle states. Interaction acts 
only when there are two or more dressed particles present. This is equivalent to saying 
that self-interaction effects are not present in the dressed particle approach. Naturally, 
there should be no dressing in our theory, because we are working with particles which 
are already fully dressed. 


Conclusions 

In this paper, we considered a model theory which provides a simplihed description 
of interactions between electrons and photons, similar to exact interactions in QED. 
We found that in order to get sensible results for the F-matrix, the theory should 
be renormalized, just as QED, by adding inhnite counterterms to the Hamiltonian. 
We also demonstrated that the renormalized version of the theory can be reformulated 
entirely in terms of dressed particles and their interactions without affecting predictions 
of the theory about scattering cross-sections or bound state energies. The dressed 
particle Hamiltonian H'^ is hnite, thus, the ultraviolet divergences are not present 
anymore. There are no bare and virtual particles in the dressed particle approach. This 
simplihes substantially the physical interpretation of the theory. From the operator 
form of the dressed particle Hamiltonian (insi) it is clear that dressed particles interact 
with each other via instantaneous forces that generally do not conserve the number of 
photons. 
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